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ABSTRACT 

We study axisymmetric perturbations of neutron star endowed with a strong magnetic 
field (magnetars), considering the coupled oscillations of the fluid core with the solid 
crust. We recover discrete oscillations based mainly in the crust and a continuum in 
the core. We also confirm the presence of "discrete Alfven modes" in the gap between 



two contiguous continua (see van Hoven & Levin (2011)) and, in addition, we can 
resolve some of them also inside the continua. Our results can explain both the lower 
and the higher observed quasi periodical oscillations (QPOs) in SGR 1806-20 and 
SGR 1900+14 and put constrains on the mass, radius and crust thickness of the two 
magnetars. 

Key words: asteroseismology - MHD - stars: neutron - stars: oscillations - stars: 
magnetic fields - gamma rays: general 



1 INTRODUCTION 

Over the last few years, a number of observational discover- 
ies have brought magnetars (ultra-magnetized isolated neu- 
tron stars) to the forefront of researchers attention. These 
extreme objects comprise the Anomalous X-ray Pulsars 
(AXPs; 10 objects) and the Soft Gamma-ray Repeaters 
(SGRs; 5 objects), which are observationally very simi- 
lar classes in many respects. They are all slowly rotat- 
ing X-ray pulsars with spin periods clustered in a nar- 
row range (P ^ 2-12 s), relatively large period derivatives 
(P - IQ-^^ - 10-^°ss-^), spin-down ages of 10^ - 10^ yr, 
and magnetic fields, as inferred from the classical magnetic 
dipole spin-down formula, of 10^ 
view see 



inate from large-scale rearrangements of the inner field or 



10''' G (for a recent re- 



Mereghetti (2008)) 



SGRs undergo periods of activity during which recur- 
rent bursts with sub-second duration and peak luminosities 
of ~ 10^^ — 10^' erg/s are observed. SGRs also show, on 
rare occasions, extreme events known as giant flares. These 
are characterized by an initial spike of duration comparable 
to that of recurrent bursts, but many orders of magnitude 
higher luminosity. Only three giant flares have so far been 
observed in over 30 yr of monitoring. 



According to the magnetar model ( [Thompson Dun- 



can|1993 2001 ) energy is fed impulsively to the neutron star 



magnetosphere when local "crustquakes" let magnetic helic- 
ity propagate outwards, giving rise to recurrent bursts with a 
large range of amplitudes. Giant flares are believed to orig- 
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catastrophic instabilities in the magnetosphere (Thompson 
Duncan|2001||Lyutikov|2003t . Most of this energy breaks 



out of the magnetosphere in a fireball of plasma expand- 
ing at relativistic speeds which results in the initial spike of 
giant flares. The decaying, oscillating tail that follows the 
spike displays many tens of cycles at the neutron star spin 
rate. This is interpreted as being due to a "trapped fireball" 
which remains anchored inside the magnetosphere and cools 
down in a few minutes. The total energy released in this tail 
is ^ lO^'^ erg in all three events detected so far. 

A power spectrum analysis of the high time resolution 
data from the 2004 Dec 27 event of SGR1806-20, observed 
with the X-Ray Timing Explorer (RXTE), led to the dis- 
covery of fast Quasi Periodic Oscillations (QPOs) in the 
X-ray flux of the decaying tail of SGR ( [Israel et al] ( |2005| )). 
QPOs with different frequencies were detected, some of 
which were active simultaneously and displayed highly sig- 
nificant QPO signals at about 18, 26, 30, 93, 150, 625 and 
1840 Hz dWatts k Strohmayer| |2006 ) . A re-analysis of the 
decaying tail data from the 1998 giant flare of another mag- 
netar, SGR 1900+14, revealed QPOs around frequencies of 
28, 54, 84 and 155 Hz (Strohmayer & Watts 2006). Hints 
for a signal at ^ 43 Hz in the March 1979 event from SGR 
0526-66 were reported as early as 1983 ( Barat et al.||1983 ). 
All QPO signals show large amplitude variations with time 
and especially with the phase of the stars rotational modu- 
lation. 

QPOs have also been argued to provide independent ev- 
idence for superstrong magnetic fields in SGRs [Vietri et ah] 
( [2007, ) . Numerous explanations have been proposed for the 
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origin of the QPOs including the torsional oscillations of 
the crust alone or even as global seismic vibration modes 



of magnetars (Piro 2005; Sotani et al. 2007a; Samuelsson 
[fcAndersson 2007 ; Steiner Watts 2009). Moreover 



Levin 



( 2006 ) argued that the QPOs may be driven by the global 



mode of the magneto-hydrodynamic (MHD) fluid core of 
the neutron star and its crust, rather than the mechan- 
ical mode of the crust. Following thi s idea, [Sotani et al.| 
(|2008); Colaiuda et al.| (2009) ; Cerda-Duran et al.| ( |20'Q9t 
recently made two-dimensional numerical simulations (both 
linear and non-linear) and found that the Alfven oscillations 
form continua which may explain the observed QPOs. The 
weak point of these simulations is the absence of the crust. 
Moreover, some of the expected eigenfrequencies (^DuncaUj 



19981 [Israel et al.|2005||Piro|2005||Sotani et al.|2007a||2008| 



Colaiuda et al.[[2009[ [Cerda-Duran et al.'"2009') match the 
observed QPO frequencies found by [Strohma yer W atts] 
(2006). These findings have opened the field of magnetar 



seismology, which provides an unprecedented probe of the 
star's crust and interior. In retrospective, predictions made 
in[Colaiuda et al.[([2009| led to re- analysis of the 2004 Dec 27 



event and revealed some new frequencies (Hambaryan et al. 



2011) 



In this work we move to a more general model for neu- 
tron star i.e. we assume a fluid interior and a solid crust 
which are threatened by a dipole magnetic field. In the in- 
terior the perturbations will be dominated by Alfven waves 
propagating along the field lines forming a continuum spec- 
trum as it has been seen in ( Sotani et al. 2008 Colaiuda et al 



2009[ [Cerda-Duran et al.[ 2009 ) while in the crust torsional 
oscillations dominate (Sotani et al.[[2007a ). Our results dif- 



fer partially from the ones by Gabler et al. (2011), where 
the authors study a magnetised neutron stars with the pres- 
ence of an elastic crust: in fact, we discover the presence of 
a type of "discrete Alfven modes" which are not present in 
the absence of a crust. These modes have been seen by [van[ 
Hoven & Levin (20TT| in the gaps between two contiguous 



continua (for this reason they call them also "gaps modes" ) . 
However, we are able to resolve these modes also inside the 
continua. 



2 DESCRIPTION OF THE PROBLEM 

Here we give a brief description of the equations and the 
boundary conditions that we used in our study. First, we 
consider a spherically symmetric and static star, described 
by the TOV equations and the line element 



ds 



dr ^ + r^(d6'^ +sin^ 6'#^) 



(1) 



Although here we used an ultra strong magnetic field 
B = 4x10^^ Gauss, we neglect the induced deforma- 
tion due to magnetic pressure, since it is still small be- 
cause the magnetic energy Em is a few orders of magnitude 
smaller than the gravitational bi nding energy Eg . Typi cally, 



Em/Eg - 10~^(B/10'^G)^ see Colaiuda et al. 



(2008) and 



Haskell et al. (2008) for results on magnetar's deformation. 
Axial and polar perturbations don't couple to each other 
and could be evolved independently, because we consider 
pure axisymmetric perturbations on a spherically symmet- 
ric star. Note that this assumption is not valid in a non- 
axisymmetric background: or example, as we mentioned be- 



fore, a strong magnetic field will slightly deform the star 
away from sphericity and in this case there is a coupling of 
axial and polar oscillations. We have chosen to ignore the 
effects of smaller deformations since although they will in- 
troduce a tiny coupling between the polar and axial modes, 
in practice the properties of the two spectra will not be af- 
fected neither qualitative or quantitative. Other effects, such 
as resistivity, can potentially have more significant contribu- 
tions (see Lander Jones (2010 )). 

In the following we consider only torsional oscillations, 
which are of axial type and do not induce density variations 
in spherical stars: in this way, the radiative part of the metric 
describing the gravitational field does not vary significantly. 
For this reason the frequency of torsional oscillations is de- 
termined with satisfactory accuracy even when the metric 
perturbations were completely ignored by setting 6g^u — 0, 
i.e. working in the Cowling approximation (Cowling ( 1941 )). 



In [Sotani et al. ( 2007a| ) the MHD oscillations of the equi- 
librium model were derived by transforming the perturbed 
linearized equation of motion and the perturbed magnetic 
induction equations into 2i^-wave equation for the displace- 
ment function 3^(t, r, 0) which is related to the contravariant 
component of the perturbed 4-velocity, 5u^ via 



(2) 



Sotani et al.[ ([2007a) for an analytical derivation . The 



2D- wave equation for the displacement 3^(t, r, 0) has the fol- 
lowing form: 



At 



Ao 



+^11^+^02 



dr^ drde dO'^ dr dO' 

(3) 

All coefficients ^tt, ^ii, ^02, ^lo and Aqi depend only 



on r and 0, see Sotani et al. (2008). In Colaiuda et al. (2009) 



a coordinate transformation of the form 

X — ^^/oisinO and Y — ^^/oicosi 



(4) 



is used to transform the 2D- wave equations into a ID- wave 
equation in the part of the star where the shear modulus 
was zero (core). The final form of this equation is: 



T d^y ~ dy 



where 
A02 - 

Aoi = 



1 /2 

aiai , 



X 



ai 



-aie 



4:7Tjie 



dV ' 



ai 
2 



(5) 

(6) 
(7) 



Here ai(r) and ji{r) are the radial components of the elec- 
tromagnetic 4-potential and the 4-current and a prime indi- 
cates derivative with respect the radial coordinate, respec- 
tively. Note that the function ai(r) is dimensionless. 

The distribution of the magnetic field inside the star can 



be found by solving the Grad-Shavranov equation ( |Grad 
Rubin[ ( [19581 , [Shafranov] ( |1966| )): 



-2A / 

e ai 



-47rji 



(8) 



with the appropriate boundary and initial conditions, see 



Colaiuda et al. (2008) for details. As discussed in Colaiuda 



et al. (2009), ai(r) shows a maximum inside the star and 
then the transformation (|4| is not any more one to one. 
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Therefore, after this maximum we choose the minus sign in 
equation Q. 

In this paper we extend our previous work incorporating 
a soUd crust in the magnetar model. The presence of the 
crust is accompanied with a non-zero shear modulus /i in 
the coefficients of the wave equation ^ . The presence of the 
shear modulus /i does not allow for dimensional reduction 
and the transformation of equation ([3| leads again into 2D- 
wave equation, in the new coordinates X and Y . The new 
equation has the following form 



, d^y , 4 I d^y 

Au ^ ^ + gYdX + 

-A dy -A dy 



dY 



where the coefficients are 
A02 



1/2 2 



47rr4 

aie^ sin^ 



cos d 



A20 
All 

Aoi 



+ M 



4 


- — sm 
V4ai 


4 


V2ai 


X 






/ di 


/2 



+ 



4 / /2 



ai 



ai 
ai 



aie 

2A 



4ai 



aie^ cos^ ^ 



sm u cos c 



(9) 

(10) 

(11) 
(12) 

(13) 
(14) 
(15) 
(16) 
(17) 
(18) 



Note that for /x = one recovers equation ([5}. The value of 
11 in the following is calculated from Vs, the shear velocity, 
see Schumaker Thorne| ( |1983| ) 



A ■ ^ 

4:7Tjie 



aie 



ai 
2 



27r Ji I cos 



4ai r 



— A cos 6* / ai 
3fI^yale — — + fi 



■ 2tiJi sin 



: COS 



Vs = {fi/py^^ 1 X 10^ cm/s. 



(19) 



2.1 Boundary Conditions 

In order to solve the equations ^ and (|9]), appropriate 
boundary conditions must be imposed. Actually, the bound- 
ary conditions in spherical coordinates were reported in 



Sotani et al. (2007a) 



Since the wave equation is given in the {X, Y) coordi- 
nates the aforementioned boundary conditions translate as 
follows: 

• regularity at the center: y{X, y) = 0; 

• no traction on the surface for the open lines: 



'dy dy 



= 0; 



• cixisymmetry 



Y 



dy_ 
dx 



a.t X = 0: 



(20) 



(21) 



• equatorial plane symmetry for i = 2 initial data i.e. 



dy 

y^ = o at y = o; 



(22) 



• equatorial plane symmetry for i = 3 initial data i.e. 

y(X,Y) = at y = 0; (23) 

• continuity of the radial function ai when the sign is 
switched in the coordinate transformation Q at ai(r) = 0. 

An additional boundary condition is needed at the base 
of the fluid-core interface: this condition in spherical coor- 
dinates is given by 

(2^- l)(2^ + 3) 



ry(-) = 



1 + 



3{& + i-l) vl 



dry^ 



+) 



(24) 



where va — B/{ATip) ' is the Alfven velocity, see 



et al.| ( |2007a| ) and |Glampedakis et al.] ( |2006t . In the {X,Y) 
coordinates the previous conditon becomes: 



Sotani 



'dy dy 
^dx^^W 



(-) 



1 + 



(2^- l)(2^ + 3) v^. 



3{P + i 



dy dy 
^dx^^W 



(25) 



(+) 



Although in this two dimensional study we have not used 
decomposition in spherical harmonics, we kept the relation 
(25) in order to be able, during the mode recycling, to infer 
in the best possible way the form of the eigenfunction that 
corresponds to a specific angular index i. For the initial runs 
we have used the more generic formula: 



1 + Vs/v\ 



r3^( + ) 



(26) 



which has been also used by Gabler et al. (2011) to extract 



the general picture of the spectrum. In general, neither the 
spectrum nor the eigenfunctions depend critically on this 
choice. 



2.2 Numerical method 



We used two representative EOSs the APR (Akmal et al. 



1998| ) and the WFF ( [Wiringa etaL| ( | 19881 ), for various mass 
models and different values of the magnetic field. At the end 
the ones that seem to fit better the observational data are 
the ones with a mass of 1.4M0 while the value for the mag- 
netic field strength which is in accordance to independent 
estimations was i?^ = 4 x 10^^ Gauss on the pole. 

We used a numerical equidistant grid 60 x 60 in 
the (X, y) coordinates, setting Xmax = \/cli max sin^ and 
^max = V^i max COS and Varying from to 7r/2. The ac- 
curacy of the code was tested, performing a simulation with 
a 90 X 90 equidistant grid: the results show that the fre- 
quencies are not significantly influenced by a change in the 
number of grid points. The stability of the scheme was also 
checked: our simulation lasted 2 seconds without showing 
any signs of instability. The base of crust (Xcrust , ^crust ) is 
taken at p = 2.4 x 10^^ g/cm^, according to the crust model 
by N egele k Vautherin] ( |1973| ) (NV). In this way the study 
is divided into two evolution problems coupled via the the 



interface condition ( 25 ) . 

First, we evolve the 1-dimensional wave equation ^ 
till the base of the crust {Xcmst, ^crust)- Special care should 
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be taken at the point where the coordinate transformation 
Q is changing sign, i.e. at ai max- In practice, one evolves 
along "magnetic strings" throughout the core of the star 
till the base of the crust. Then in the crust, we evolve the 
2-dimensonal wave equation (|9|, and the oscillating "mag- 
netic strings" of the core are replaced by a "magneto-elastic 
membrane" . 

As a test run for our model, we reproduce the results in 



Sotani et al. (2007a) for crust oscillations. In this paper, the 



authors study the torsional oscillations of a non-magnetized 
star. They also consider, as a first approximation, the no- 
traction condition in the core-crust interface, i.e. instead of 



the condition ( 24 ) they require that 3^ has to be continuous 
through the interface i.e. 3^^^^ = -^(-)- H^^^? we cannot set 
the magnetic field equal to zero since our coordinate system 
(X, Y) is based on the magnetic field strength i.e. on ai(r). 
However, it was shown in Sotani et al. (2007a| that the in- 
fluence of the magnetic field on the frequencies becomes im- 
portant only if B > 10^^ Gauss, thus we can choose a very 
low magnetic field for magnetars, e.g. B = 10^^ Gauss, and 
perform test runs in order to reproduce the results of'SotaniJ 



et al. (2007a). Using this magnetic field and the no-traction 



condition on the interface, we find the results reported in 
Sotani et al. (2007a). This test verified that the code can 



reproduce earlier results for the torsional oscillations of the 
crust. 

In a similar way we can isolate the continua found in 
( 2009| and trace its changes due to the pres- 



Colaiuda et al. 



ence of the crust. One of the characteristic properties of the 
continua is the scaling of the frequencies that seem to form 
the edges of the continua. That is, the edges of the various 
continua follow the rule that has been found for odd type of 



initial data in Colaiuda et al. (2009) i.e. that 



fLr, ^ {n+ 1)/lo and fur, ^ (n + l)fuo, 



(27) 



where /l stands for the lower frequency of the continua with 
an eigenf unction located mainly near the polar axis, and fu 
stands for the upper frequency of the continua with maxi- 
mum amplitude near the last open magnetic field line. These 
scaling laws for the continua can be easily spotted in Figs. 
□ and [a 

In this coupled core-crust problem apart from the con- 
tinua and the crust modes, we found the imprints of a new 
family of discrete modes, a detailed description of these 
modes is presented in the following section. 

Finally, we made some additional numerical tests i.e 
we set initial data in order to excite only the crust and to 
study the propagation of the energy in the system. We pick 
up two points, one in the core and one in the crust, and 
we follow the evolution of the perturbation during the time 
of our simulation, performing a FFT at different times and 
comparing the amplitude of their peaks. We found that the 
energy quickly flows from the crust towards the core, ex- 



citing the Alfven continua as seen by Levin (2007) and by 



Gabler et al. (2011). In the same way, the perturbations 



of the core excite the crust: for this reason it is more ap- 
propriate to define the modes generated by the continua, 
global modes, i.e. modes that involve both crust and core 
oscillations. The conclusion of this test is in agreement with 
earlier suggestions ( Glampedakis et ar]|2006 ) claiming that 
when the magnetic field permeates the whole the star, the 
perturbations cannot be confined in the crust. Another no- 
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Figure 1. Identification of the frequencies of SGR 1806-20. We 
show that a stellar model APR with mass M = IAMq and radius 
R = 11.57km can explain all the frequencies observed. Upper 
panel: the magnetic field strength is 4 x IQ-'^^G. Lower panel: the 
magnetic field strength is 2 x IQ-'^^G. In both panels the blue bands 
represent the continuous spectrum of the core oscillations. 



table result is that even when we associated a mode to the 
crust, the actual eigenfunction was not confined in the crust 
but it was extended throughout the star, see Fig. [3] In ad- 
dition, we find that when a crustal mode is embedded in 
the continua, its eigenfunction shows the same structure of 
the continua that hosts it. Contrary, when a crust frequency 
is found in the gap between the continua, its eigenfunction 
shows a structure more similar to the one of the fundamen- 
tal crustal frequencies, compare the upper panel of Fig. |3] 
and the lower panel of Fig. [4] 



3 RESULTS & DISCUSSION 

The scope of this paper is not only to understand if the 
oscillations are global modes (i.e. combined crust-core os- 
cillations), pure crustal modes or pure core modes but also 
to identify the magnetar models that can fit the observed 
QPOs. This can potentially lead to a better understanding 
of the order with which the various modes were excited and 
left their imprints in the signal. Moreover, a proper identi- 
fication of the observed QPOs will unveil the details of the 
magnetar structure. 
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crust APR14 B = 4.25 X 10^^ Gauss 

discrete Alfven modes 
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Figure 2. Identification of the frequencies of SGR 1900+14. Since 
just few frequencies were identified for this SGR, we show that a 
stellar model APR {upper panel) with mass M = 1.4M0, radius 
R = 11.57km and magnetic field B = 4.25 x lO-'^^G can explain 
all the frequencies observed, as well as a stellar model WFF with 
mass M = 1.4M0, radius R = 10.91km and magnetic field B 
1 X lO-'^^G {lower panel). In both panels the blue bands represent 
the continuous spectrum of the core oscillations. 



3.1 Analysis of the spectrum 

As we already pointed out and as it is shown in Fig. [l] and 
Fig. |2] the spectrum of the oscillations is composed by three 
different types of modes: 

crustal modes: they are associated with the crust and 
they acquire a discrete spectrum. We found that, when the 
magnetic field increases, some modes are absorbed by the 
continua, compare upper panel of Fig. [H] where B ^ 10^^ 
Gauss, with the lower panel of the same figure, where 
B c:^ 6 X 10^^ Gauss: it is clear that some modes disap- 
pear. The strong coupling with the crust leads to an effec- 
tive transfer of energy from the crust to the core, especially 
when the crust frequency is embedded in the continua (see 
for example the 66Hz frequency in the upper panel of Fig. 
[T]): from a combined analysis of the peak amplitude of a 
crust frequency embedded in the continua and of the peak 
amplitude of the edge of the hosting continua, we found that 
the energy lost by the crust frequency is stored in the edges 
of the continua. In addition, we found that the fundamental 
crustal frequency does not scale with the increase or de- 




Figure 3. Typical eigenfunctions of a "crust" mode for the model 
APRi4 that identifies the SGR 1806-20 (upper panel) and the 
model WFFi4 (lower panel), that identifies the SGR 1900+14. 
This latter frequency is unlike to be observed, since it does not 
reach the stellar surface. The noise visible in both figures corre- 
sponds to the change of sign of the (X, Y) coordinate at the point 



crease of the magnetic field: this behaviour assures us that 
the crustal frequencies are not scaling with the magnetic 
field unlike the Alfven modes. For magnetic field B > 10^^ 
Gauss, the relation among the fundamental frequency and 
the overtones, is given by: 



h 



\ />crust 



for l> 2 and 



0. 



(28) 



For a magnetic field B < 4 x 10^^ Gauss, the law found by 
Samuelsson & Andersson] ( |2007 |) is recovered: 



h 



V(^-l)(^ + 2) 



2tiR 



for 



0. 



(29) 



For intermediate magnetic fields, a proper scaling cannot 
be found. The problem of the scaling of the crustal mode 



was treated also by Gabler et al. (2011), where the authors 



discuss the appearance or disappearance of the crustal mode 
in relation to the magnetic field strength. 
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Figure 5. Upper panel: the crust al frequencies of an APR14 
model for a weak magnetic field B ~ lO-*^^ Gauss. Lower panel: for 
a weak magnetic field 5 ~ 6 x lO-*^^ Gauss: some of the frequencies 
disappear as the magnetic field increases. 



Figure 4. Eigenfunctions of the "crust" mode 66Hz for the 
APR14. Upper panel: the the magnetic field strength is B = 
4 X 10^^ G and the frequency is hosted in the continua (model 
for SGR 1806-20). Lower panel: the frequency is located in a gap 
and the magnetic field is B = 4.25 x lO^^G (model for SGR 
1900+14). 



continuous spectrum: the continuous spectrum is gen- 
erated by the presence of a strong magnetic field (B > 10^^ 
Gauss). It is generated in the magnetised fluid present in the 
core but, in the case of strong magnetic field it extends and, 
through the coupling with the crust, its oscillations pene- 
trate also in the crust, giving rise to global modes. When 
the magnetic field increases, the continuous spectrum be- 
comes more and more energetic, absorbing also some crust 
modes. In addition, the edges of the continua scale with the 
magnetic field, since / = B/ (^47rpL), where L = L(X) is 
the length of the "magnetic strings" along which the per- 
turbation propagates (compare the upper and lower panel 
of Fig. [1]). The frequencies of the continua show the scaling 
described by equation (27). Note that the continuous spec- 



trum always shows gaps in its structure for the model APR 
(see Fig. [TJ, while the model WFF does not have a gap be- 
tween the fourth and fifth continua that overlap (see Fig. 



[2]) : the presence/absence of gaps seems then to be a char- 
acteristic of a particular model more than a general one, as 
pointed out in |van Hoven & Levin| (|2011|. 

discrete Alfven modes: these modes, as the crustal 
modes, are discrete but we identify them as Alfven modes, 
because they scale with the magnetic field as the frequencies 
of the continua. Their structure is similar to the one of the 
crustal mode (compare the upper panel of Fig. [t] and Fig.js] 
with Fig. [3]) but their scaling follows the scaling observed in 
the continua: 

/i°)«(n + l)/o(°). (30) 
Note that those modes have not been observed in the case of 



the absence of a crust (see Colaiuda et al. (2009)) . For this 
reason, they can be interpreted as an effect of the coupling 



between the fluid core and a solid crust (see van Hoven & 
|Levin| ( |2QTT| ). 

3.2 Identification of the QPOs 

From the variety of magnetars models that we have exam- 
ined in order to fit to the observational data only a few pro- 
vide oscillation frequencies that are in agreement with the 
observed QPOs. More specifically, we compared the data 
that our numerical code produced with the frequencies from 
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Figure 6. Eigenfunctions of the first two edges of the first con- 
tinuum for the model APR14 (SGR 1806-20). 



the timin g analysis of the SGR 1 806-2 and SGR 1900+14 
shown in Strohmayer & Watts (2006). In this paper, 



the 

authors identify several QPOs, of different duration, in the 
tail of the two events. In particular, for the more recent 
event, SGR 1806-20, the identified frequencies are 18, 26, 
30, 92, 150, 625 and 1840 Hz while for the SGR 1900+14 
they found the following frequencies: 28, 53, 84 and 155 Hz. 



A more recent study by|H ambaryan et al. (2011), based on 
predictions by Colaiuda et al.^ ( ^2009 ), confirms, by using a 
different analysis technique, the earlier results and, in addi- 
tion, unveils the presence of at least three new frequencies 
for the SGR 1806-20 these are 16, 21 and 36 Hz. 

The new results pose extra challenge, since already the 
previous studies of pure torsional oscillations in magnetars 
could not explain all observed frequencies, in particular the 
lower ones (18, 26, 29 Hz), mainly because of the small spac- 
ing among them. If the two new frequencies (16Hz and 21Hz) 
are added then it is impossible by any means the expla- 
nation via crust oscillations alone. On the other hand, the 




Figure 7. Typical eigenfunctions of a "discrete Alfven mode " for 
the model APR14 (SGR 1806-20). Note that the 53Hz frequency, 
represented in the lower panel, is unlike to be observed, since it 
does not reach the stellar surface. 



but still both calculations did not take into account the pres- 
ence of a crust. 

Here we show that our modelling of magnetar dynam- 
ics can explain all the observed frequencies, and inversely 
via this explanation we can constrain the parameters of the 
magnetar (radius, mass, equation of state and magnetic field 
strength). In Figs. [l] and [2] we graphically show our findings, 
which comprised by 3 types of oscillations, the Alfven con- 
tinua in the core, the discrete Alfven modes (whose eigen- 
functions are plotted in Figs, [t] and |8| and the crust al tor- 
sional modes confined mainly in the crust (whose eigenfunc- 
tions are plotted in Fig. |3|. The eigenfunctions were found 
using the eigenfunctions recycling program already used in 
Gaertig k Kokkotas] (|2009) (see there for more details). 

A model with the equation of state (EOS) APR for the 
core and the NV model for the crust, with mass M — 1.4M0 
and radius R = 11.57km, could fit the new data discovered 



Alfven continua in the core (Colaiuda et al. 2009 Cerda- 



Duran et al. 2009) offer better chances for an explanation 



by Hambaryan et al. (2011 ) as well as the previous data for 
the SGR 1806-20 ( [Watts Strohmayer|[2^ . Here, we fo- 
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Figure 8. Typical eigenfunctions of a "discrete Alfven mode 
for the model WFF14 (SGR 1900+14). 



cus on the first five frequencies, i.e 16, 18, 22, 26, 29 Hz, since 
they are the more difficult to be explained, the higher ones 
can be explained in various ways, as multiples of the lower 
ones or as polar modes (see Sotani Kokkotas (2009)). 

In the upper panel of Fig. [6j the eigenfunction of the 
16Hz frequency, i.e. the first edge of the continua (see Fig. 
[1]), is shown. The perturbation involves both core and crust 
and it can be classified as a global oscillation, with the crust 
and the core oscillating together due to the strong coupling 
induced by the magnetic field. The perturbation seems to be 
localised near the magnetic axis and just outside the region 
with the closed lines. In the bottom panel of Fig. [6] the 
eigenfunction of the second edge of the first continuum, the 
18Hz frequency, is shown: its structure is similar to the one 
of the 16Hz frequency, but it shows broader spreading and 
intense oscillation. 

The 22 Hz mode has a discrete nature, as we already 
pointed out in the previous section, see upper panel in Fig. 
[3] Studying the time evolution of the perturbation, as ex- 
plained in the previous section, we find that the core is ex- 
cited after some time, when the crust successfully transfers 



its energy to it. It seems to be a fossil of crust mode found 
for weaker magnetic field and it is then reasonable to clas- 
sify it as crust al frequency. The next frequency found, at 26 
Hz, see Fig.jT] is similar in structure to the crustal one but, 
analysing its time evolution, we find that it originates at the 
crust core interface and its nature is discrete. In addition, 
although it is a discrete frequency, it scales with the mag- 
netic field: it has then all the characteristics of a discrete 
Alfven mode. 

Finally, at a frequency of ^ 30 Hz, we observe excita- 
tion mainly of the core region but we notice also significant 
excitation in the crust and thus this frequency can be iden- 
tified as a global oscillation, that is localised mainly in the 
core and then forces the crust to oscillate violently. 

Note that the higher frequencies are located closer to 
the magnetic axis than the lower ones. This behaviour was 
already seen in Colaiuda et al. ( 2009) and Cerda-Duran et al.| 



(2009). Note also that the part of the star that seems not 



to be excited corresponds to the closed magnetic field lines. 
As we already noted in Colaiuda et al. (2009), the closed 



magnetic lines have a significant smaller amplitude than the 
open magnetic field lines and cannot be observed via oscilla- 
tions or thermal surface phenomena, since they are confined 
in the interior of the star. 

Concerning the SGR 1900+14, an unique model that 
can fit better the observed frequencies cannot be found be- 
cause of the paucity of the observations. For this reason, 
we find that two models can fit the identified QPOs for the 
SGR 1900+14: the EOS APR with a mass M = 1.4M0 
, radius R = 11.57km and a magnetic field strength B = 
4.25 X lO^^G, as well as the EOS WFF with mass and ra- 
dius respectively M = 1.4 M© and R = 10.91km, and a 
magnetic field strength B = 4 x lO^^G. We focus our study 
on the identification of the first three frequencies observed: 
28Hz, 53Hz and 84Hz. In the case of APRi4with a magnetic 
field strength B = 4.25 x lO^^G, all the three frequencies 
can be identified as discrete Alfven modes, see upper panel 
of Fig. [2] The oscillations involve both crust and core. Con- 
trary, in the case of model WFF14 (see lower panel of Fig. 
[2]), the frequencies can be identified as global modes (the 
28Hz and the 54Hz) and as crustal mode (the 84Hz). In 
particular, the 28Hz frequency is identified as the first edge 
of the second continuum while the 54Hz is identified as the 
second edge of the third continuum (see Fig. [2]): both those 
frequencies are global modes, with excitations that involve 
both the crust and the core. The 84Hz is a crust frequency: 
its nature is discrete and the oscillations, although initially 
confined in the crust, expand rapidly in the core. Also in 
this model, we find a discrete Alfven mode around 26Hz. 

Note that the equation of state used for the crust in 
all the models that we present in this paper is the NV. As 
already pointed out in 'Sotani et al. ( 2007b|), the use of the 
EOS proposed by Douchin & Haensel| (|2001|) for the crust 
does not alter significantly the results. Still since it produces 
a thinner crust, it creates difficulties in the simulations since 
it demands finer grid in the crust in order to evolve properly 
the perturbations inside it. 

The results present in this section partially agree with 
the ones found by Lee ( 2008 ) , where a similar problem of a 



star with solid crust and fluid core both permeated by a mag- 
netic field is studied in Newtonian theory as an eigenvalue 
problem. Without time evolutions only discrete Alfven fre- 
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quencies could be found and thus no continuous spectrum 
appears. However, the discrete Alfven frequencies that we 
found here agree quahtatively with the ones described by 
Lee. 

It is worth mentioning that we analysed also alterna- 
tive combination of magnetic field configuration, consisting 
of both poloidal and toroidal magnetic fields. We found that 
the presence of a toroidal component in the background 
magnetic field shifts the spectrum towards lower values. This 
extra parameter complicates significantly the procedure of 
identifying the magnetar model that will better fit the ob- 
served frequencies. 

Note that a more accurate treatment of the shear mod- 



ulus than the one used in equation (19) may change quan- 
titative the crustal modes. In addition, the presence of a 
superfluid component in the inner crust and/or in the core 
could affect the frequencies significantly, as discussed in|An-| 



dersson et al. (2009). 



4 CONCLUSIONS 

We studied the torsional oscillations of a magnetar in a gen- 
eral relativist ic framework, consisting of a relativist ic neu- 
tron star with a solid crust and a fluid core. The core and 
the crust are coupled by the strong magnetic field and this 
coupling is defined by the appropriate boundary conditions. 
We find that the presence of a crust and its coupling with 
the core partially alters the earlier picture presented in |Co-| 



laiuda et al. (2009). 



In particular, the presence of a crust makes the spec- 
trum denser and thus offers an explanation of the nature of 
all the observed frequencies. We can distinguish, using the 
eigenfunction, between global modes (i.e. modes for which 
both the crust and the core oscillate), crust modes and dis- 
crete Alfven modes (i.e. modes that have a discrete nature 
but that scale with the magnetic field) . Unlike the paper by 
van Hoven & Levin ( |2011 ), we find that also in the case of 
intermediate magnetic field {B < lO^^G), the oscillations is 
not confined to the core but also the crust is excited. We find 
a family of modes, the discrete Alfven modes, in the gaps 



between two contiguous continua (as found by van Hoven 
& Levin (20111)) and, in addition, we can resolve them also 



inside the continua. 

We can identify all the frequencies observed in SGR 
1806-20 and SGR 1900+14 and this allows us to put some 
constrains on radius, mass, crust thickness and magnetic 
field strength. In particular, the frequencies observed in SGR 
1806-20 can be fitted uniquely with the APR stellar model 
with a mass M — lA M©, a radius R — 11.57km, a com- 
pactness M/R = 0.178 and a crust thickness Ar/i? = 0.099. 
The magnetic field strength at the pole is B — A x lO^^G. 
Contrary, the SGR 1900+14 cannot be strongly constrained 
because of the limited information from the observations, 
i.e just three QPOs were identified. For this reason, the 
SGR 1900+14 cannot be reproduced uniquely by a single 
model, as it has been done for the SGR 1806-20, since at 
least two stellar models could reproduce its observed fre- 
quencies with a good accuracy: the APR stellar model with 
a mass M = 1.4 M©, a radius R — 11.57km, a compactness 
M/R = 0.178 and a crust thickness Ar/i? = 0.099 but with 
a magnetic field strength B = 4.25 x lO^^G and the WFF 



model with a mass M — 1.4M0, a radius R — 10.91km, a 
compactness M/R = 0.189, crust thickness Ar/R = 0.085 
and a magnetic field strength B = 4 x lO^^G. Note that, in 
all the models, the crust is described by the NV equation of 
state. 
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